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Abstract 

We analyze here late evolutionary stages of massive (Mo > 8Mq ) close binary 
stars. Our purposes are to study possible mechanisms of gamma ray bursts (GRBs) 
origin. We suppose in this paper that GRB phenomenon require formation of mas- 
sive (~ IMq) compact (R < 10 km) accretion disks around Kerr black holes 
and neutron stars. Such Kerr black holes are products of collapse of Wolf-Rayet 
stars in extremely close binaries and merging of neutron stars with black holes 
and neutron stars with neutron stars in close binary systems. Required accretion 
disks also can be formed around neutron stars which were formed during collapse 
of accreting oxygen-neon white dwarfs. We have estimated frequencies of events 
which lead to a rotational collapse concerned with formation of rapidly rotating 
relativistic objects in the Galaxy. We made our calculations using the "Scenario 
Machine". 



1 Introduction 

Investigation of the gamma ray bursts (GRBs) physics is one of the most actual as- 
trophysical problems during last decades and it's popularity permanently grow. Some 
years ago frequency of appearance of the articles on GRB physics had exceeded the fre- 
quency of the recorded GRBs which equals to approximately one flash per day for the 
detectors with threshold level ~ 1CP 7 erg cm -2 sec -1 in the range 30 KeV - 100 MeV 
1641 1721 . About thirty year have elapsed since pioneering works about observed GRBs 
I37ll96ll59l . but reliable observational limits on the main parameters of the events have 
appeared only during recent years and two conceptions which probably present two 
main types of GRBs were accepted. These types are short (SGRBs) which duration is 
lower than 2 seconds and long (LGRBs) which duration is ~ 2 — 200 seconds | 38ll34l . 
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GRBs were supposed to have cosmological nature by Usov and Chibisov 1 96 1 . Now 
numerous GRBs positions were identified with positions of galaxies and some super- 
novae type I b,c since 1997 1 18 1. So, cosmological nature of the most of the GRBs had 
become evident |96 65 66|[M][72|[21|34|. It can be illustrated by the isotropy of the 
GRBs distribution in the sky and by the deficiency of faint flashes in comparison with 
isotropic distribution in the Euclidean space [51 82, 83l ll5l . These facts enable re- 
searchers to realize cosmological distances to the most GRBs and to reconstruct some 
cosmological parameters, also they put problem of reconstruction of the star formation 
history in the Universe |89 15 1. Direct identification of red shifts of some GRBs (until 
z ~ 6) [23 1 supports this possibility. Identification of the host galaxies enables to es- 
timate distances to some GRBs and their total energy in the assumption of spherically 
symmetric radiation - 0.1M Q c 2 172 1. and frequency of the GRBs in a galaxy with mass 
like the mass of the Milky Way ~ 10~ 6 - 10~ 7 yr -1 (58||67|. Assumption about col- 
limation of the gamma radiation in the narrow space angle (^0.1 steradian) 1671 1571 
allows to get the last estimations up to the energy ~ 10 51 — 10 52 erg and frequency 
~ 10~ 4 — 10~ 5 yr -1 . Note that these numbers are still uncertain by a factor of ~ 10. 
Also it is important to note that GRBs frequency is less than frequency of the known 
supernovae at least several tens times. This circumstance directed investigators to the 
most relic, but powerful events in stellar life. It is important to say that energetics of 
SGRBs in X-ray and gamma ray ranges is almost one hundred times less than ener- 
getic of LGRBs 1171 . This is direct indication that frequency of SGRBs per galaxy 
can be higher than frequency of LGRBs. Additional indication of relatively higher 
frequency of SGRBs is relatively small distance to typical identified SGRB (z « 0.2) 
with comparison to the distance to typical GRB (z=2) |42|). 

Millisecond variability of observed flux of GRBs is evidence of a small volume 
of the main energy-release region, it's size does not exceed ~ 10 8 cm. Three sorts 
of astrophysical objects with such dimensions are known: degenerated dwarfs (DD, 
or white dwarfs, WD), neutron stars (NS) and black holes (BH). High energy release 
~ 10 51 — ~ 10 52 erg is typical for NSs and BHs. It provides for conditions upon ideas 
how to explain GRB mechanism. Concept of merge of two NSs under influence of 
gravitational waves was advanced to depict SGRBs. Roche lobe overflow by one of 
the compact stars leads to it's decay in dynamical time scale 10~ 4 — 10~ 3 sec. |91 1. It 
allows to conclude that energy release during such event is enough to produce a SGRB 
and to make a relation between a SGRB and a NS+NS merging |6|. NS also can be a 
component in a close binary including BH. List of the four known close binary systems 
consisting of neutron stars and radio-pulsars in the Galaxy which are able to merge in 
timescale shorter than Hubble time was compiled |32|. It includes such systems as 
B1913+16 with orbital period about eight hours. According to the scenario program 
such pairs also have to merge with frequency equals to frequency of merging of the 
close binary systems 1 • 10~ 4 yr -1 1 501 1921 19~3*1 . In the present article we consider 
both variants of SGRB formation to estimate their frequencies. It is necessary to note 
that in spite of formal high "delay" of merging after the moment of the formation of 
the system (NS+NS, NS+BH) considerable part of merging happen in first 1-2 Gyr. 
Although short GRBs can be found also in old elliptical galaxies 1701 1221 . most of 
them have to be related to galaxies with active current star formation |25 1. In fact one 
short GRB was found in the galaxy with active star formation |74|. Observation of 
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the short GRB 050709 in optics allows to exactly exclude even a faint SN Ic in the 
same place \27\. But it is important that in elliptical galaxies there are only SGRBs 
E3I1I3EI1IH1IZ1- It is evidence in favour of the model of NS+NS or NS+BH mergers. 

LGRBs are concerned with collapses of rapidly rotating nuclei of massive pre- 
supernova which produce Kerr BH |67 1. Collapse of a fast rotating star as mechanism 
of supernova explosion was suggested by |5|. Numerical gas-dynamic model of such 
event was constructed by |95|. There are two possible causes of fast pre-supernova 
nucleus rotation: acceleration of the star nucleus rotation by it's grip with angular 
momentum conservation 1 87 56 1 or companion presence near pre-supernova helium 
Wolf-Rayet star ( WR) which is one of components in a close binary 1301 11031 l67l 1731 
89 90 1 . The last version seems to be preferable due to very probable significant slowing 
down of rotation of cores of single massive stars during their evolution 1 88 , 60 1 . Orbital 
period of the system consisting of pre-supernova type I b,c (WR star) and another 
component has to be shorter than 1-3 days to produce Kerr black hole. At least one 
binary including WR-star (progenitor of GRB) is known at present time: Cyg X-3 1 89 1 . 
There are also some known binary systems including post-Kerr BH, for instance, V 5 18 
Per II89I . Observational data allow to make relation between LGRBs and explosions 
of SN I b,c which mean the end of WR-star evolution 1 12 16 1 in galaxies with active 
current star formation |22 A companion of a WR-star (a progenitor of a Kerr 

BH) can be a main sequence star, a BH or a NS. We investigate here such systems to 
estimate frequency of formation of Kerr black holes in them. Note that we treat a black 
hole as Kerr BH in this paper if Kerr parameter 

where I is the moment of inertia of the black hole, ft is its angular velocity, Mbh is 
the mass of the black hole. 

It is known that even in very close binaries with helium non-degenerated progen- 
itors which mass is about 2.5M NS formation does not allow to produce NS with 
period of rotation shorter than ~ 0.05 sec. according to elementary estimations 1 . At 
the same time the ultimate lower limit of the period of a neutron star is much shorter 
and equals to ~ 0.001 sec. It excludes such NSs from list of progenitors of the GRBs. 
Nevertheless probably there is another channel of formation of a fast rotating single 
NS due to merging of close binary degenerated dwarfs. And at least one of the com- 
ponents in the case has to be oxygen-neon (ONe) white dwarf. Oxygen burning in it 
does not produce enough energy to destroy compact dwarf and finally leads to collapse 
and formation of a neutron star I6l1 1631 . Initial mass of a star in a close binary sys- 
tem has to be ~ 8 — 10Af© to form oxygen-neon WD 1 30 1 . Collapse of such dwarf 
during merging of the components of the close pair will guarantee NS formation with 
over critical rotation that leads to formation of a LGRB according to our model. In 

'For Mjj e = 2.5Mq practically filling its Roche lobe with radius R[fe = 0.340 1941 major semi-axis 
a ~ i?Q for another companion of solar mass. P or j, = 10 4 M °'_^ M2 sec - ( a = o,/Rq, M = M/Mq), 
radius of the neutron star is -Rjvs = 10 6 cm., radius of the iron core of the helium star is Rp e = 3 ■ 10* 
cm H|. P NS = 10*^ (^f = 10^ ( ^) 2 « 0.06 sec. 
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general this model is similar to the model of collapse of a fast rotating WR-star and it 
is probably possible to use this model to explain "long" flashes of gamma-rays. 

At the same time compact accretion disk can be formed near young NS and super- 
critical accretion onto this NS can leads to a formation of a transient source of high 
energy photons [90|. A close companion of ONe WD has to be helium, carbon-oxygen 
or ONe dwarf. And although observational evidences of such merge stay not totally 
clear we have included these events in our analysis of frequencies. There are some 
known observable analogues of close binary degenerated dwarfs which are necessary 
to realize such scenario. They were found in the last years during search for progenitors 
of supernova stars type la. Most of known close binary degenerated dwarfs have com- 
mon mass lower Chandrasekhar limit, but three of them have masses higher than this 
limit and only one of them can merge during period of time ~ 10 10 years [62]. This 
circumstance and scenario estimations give us basis to hope that it is possible to pro- 
duce binary systems consisting of ONe WD and degenerated companion with common 
mass higher than Chandrasekhar limit which also are able to merge under gravitational 
waves influence during Hubble time. This is third possible scenario explaining obser- 
vation of "long" flashes of gamma rays. 

It is necessary to note that although the fact of possible formation of NS during ONe 
WD collapse was mentioned 1 77 1 more than once, but the picture of the explosion of 
accreting matter ONe white dwarf with NS formation is not known in details yet. Ob- 
servational manifestations of degenerated ONe dwarfs with comparable masses have 
not been detailed analyzed in spite of relatively high possible frequency of such events. 
According to scenario model it is ~ 0.01 per year in the Galaxy |93 1. But the most part 
of them are helium WDs. Duration of the phase of the destruction of a dwarf is only 
some seconds. After NS and massive disk formation the disk will evolve in dissipative 
time scale depending on the disk thickness and this time will be 10 2 - 10 4 times longer 
than Keplerian time of the disk |79 90 1. That is, process of destruction of degenerated 
companion can manifest as powerful X-ray flash with duration of some hours or days. 
Non-recurrent events of such type are actually recorded |1 78 1. It is evidently that 
observed X-ray bursts present a heterogeneous group 1 80 1 and it is a problem for the 
distant future to pick out among them events of dynamical disruption of degenerated 
components in binary stars. It is important for investigations of GRB models to dis- 
cuss how we can explain production of short flashes with magnitude comparable with 
SN explosions in close binary systems. Current classification of GRBs based on their 
duration is probably not exhaustive. Investigations of spectra of GRBs allow, for ex- 
ample, to outline indications of third subgroup of the phenomenon; it has duration 2-10 
seconds and relatively soft gamma-spectra. It is obtained that optical afterglows during 
first ~ 100 days have bimodal behaviour 1461 . This result is based on investigation of 
about forty light curves of optical "echoes" of bright GRBs |46|. Subgroup of short 
GRBs is also heterogeneous |84|. It is possible that following detailed study of these 
phenomena will allow to determine new sorts of GRBs which will differ in models of 
progenitors, in parameters of the models or jet axis orientation relative to observer. 

Formation of young neutron star due to ONe dwarf collapse with mass higher than 
Chandrasekhar limit seems to be interesting not only as potential mechanism of GRB 
formation, but also as possible way of magnetars - X-ray pulsars formation |49 1. Only 
one of them - XTE J18 10-197 was revealed as transient radio-pulsar. Magnetars are 
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neutron stars with magnetic fields ~ 10 15 Gs and rotational periods 5-12 sec. and 
current estimation of their birth frequency ~ 10~ 3 per year ll85l l39l l86l 1751 . It is 
well known that magnetic fields of some degenerated dwarfs reach to ~ 10 8 -10 9 
Gs |99|. During collapse of such dwarf with radius ~ 10 9 cm into NS with radius 
10 6 cm trapped magnetic field will be intensified by a factor of million, i.e. will 
reach to ~ 10 14 -10 15 Gs (observed value, see |102| for details). Slow rotation of 
magnetars is probably consequence of their very strong field which rapidly increases 
their rotational periods from initial, which are probably very short, up to observed 
value 0. To accelerate ONe dwarf rotation to the quantity enough to form Keplerian 
disk around young neutron star after dwarf collapse, i.e. GRB, it is enough accord- 
ing to simple estimation to accrete a part of matter of the Keplerian disk with mass 
AM/Mq ~ (Rns/Ronc)^ ~ O.O3M , where R NS and R Ne - radii of NS and 
degenerated ONe dwarf correspondingly. 

2 Population synthesis 

We use the "Scenario Machine" to estimate frequencies of the described above events 
which can lead to GRB formation. For every set of the initial parameters we have 
conducted population synthesis of 10 6 binary systems. 

Since the "Scenario machine" working principles were described more than once, 
in present work we will limit ourself only with mentioning of the most important pa- 
rameters which influence on the results of numerical modelling binaries under inves- 
tigation. Detailed description of the "Scenario machine" may be found in the book 
written by |53|. 

2.1 Initial mass ratio distribution 

Population synthesis was conducted for two types of the initial distribution of the mass 
ratio of the components in a binary system f(q) = q aq : flat a q = (the most probable 
value, see 1 40 1 for details) and quadratic a q — 2 (see 1 52 1 for details). As the coefficient 
of the mass ratio in a binary system q = m-z/mi we assume the ratio of the mass of 
the secondary companion mi to the mass of the primary companion m\ of the system 
(mi > m-2). Note that for low mass systems (mi < 10A/ Q ) we take a q — in all 
scenarios. 

2.2 'Kick' during supernova explosion 

It was supposed in our calculations that neutron star or black hole during supernova ex- 
plosion is able to get additional "kick" v, it's velocity distributed by Maxwell function: 

2 v 2 

L e -^ (2) 

it's direction is equiprobable. But the quantity of dispersion vq of the remnant is one 
of the crucial parameters for estimations of the frequencies. We should to say that the 
results of population synthesis are highly dependent on the quantity of the parameter 



5 



vq. Increasing vq higher than orbital velocity in close binary systems ~ 100 km s 
leads to sharp reduction of number of systems including relativistic companion. 

Let us suppose that absolute value of kick velocity during BH formation depends 
on the part of the mass loss by an exploding star during supernova explosion. In our 
calculations we assumed that during supernova explosion a star lose a half of it's mass 
(see |8| for details). If we assumed characteristic velocity of the neutron star kick vq, 
the quantity of the parameter vq for the case of formation of the black holes is defined 
as VQ h = 0.5vq in present work. 



2.3 Stellar wind 

Mass loss by an optical star during its evolution is still poorly known at present time. 
Although it is possible to significantly reduce uncertainties (see, for instance, JSJ), 
there are no ample reason to choose single scenario of mass outflow by stellar wind as 
a standard, so population synthesis was done using two different scenarios of mass loss 
rate by a non-degenerated star. Let us call them A and C. 

Stellar wind significantly influences on evolution of the massive stars which core's 
collapse can lead to a GRB formation. In our conception the phenomenon of GRB can 
happen only in a close binary system. It's components either merge with Kerr black 
hole formation or rotational collapse of WR-star (into Kerr black hole) or white dwarfs 
(into neutron star) is enabled due to orbital motion. Optical star mass loss greatly 
influences on the major semi-axis of binary system and is able to significantly change 
number of close binary systems which can produce GRB. 

In this work we use quasi-conservative mass transfer (see 1531 1981 for details). In 
this case we calculate the major semi-axis of the system using formula 




(3) 

where a/ is the final major semi-axis, a% is the initial major semi-axis. In this equation 
qf and <ft are final and initial values of q = M accr /Md onor , here M accr is the mass of 
the accreting star and Mdonor is the mass of the donor star. Parameter /3 = — (M accr — 
Ml^) I (M r donor — Mdonor) we calculate in this case as minimal value between (3 = 
1 and the ratio f3 — Tkh {donor) /TKH{iccr), where Tkh (donor) is the Kelvin- 
Helmholtz time for the donor, Txiiiaccr) is the same for the accretor. 

Scenario A has a weak stellar wind. The mass loss rate M during the main sequence 
(MS) stage (33 1 is 

M ~ L/Voo, (4) 

where L is the luminosity of the star and Voo is the wind velocity at infinity. 
For giants we take a maximum between (|4j and the result obtained by 1 44 1 : 

M ~ L 1A2 R a 61 /M°", (5) 
where R is the stellar radius, M is its mass. 
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For red supergiants we take a maximum between @ and Reimers's formula |41 1: 

M ~ LR/M, (6) 

The mass change AM in wind type A during one stage (except WR-stars) is no 
more than 0.1(M — M core ), where M is the mass of the star at the beginning of a 
stage and M core is its core mass. Mass loss during the Wolf-Rayet (WR) star stage is 
parametrized as 0.3 • Mwr, where Mwr is the maximum star mass during this stage. 
For calculations of stellar wind type A we used the core masses obtained by [ 1 00 1 and 

MM- 

In scenario C the stellar evolution model is based on the results of 1971 . which 
reproduce most accurately the observed galactic WR star distributions and stellar wind 
mass loss in massive stars. Calculations of mass loss by a star were conducted if we 
used the formula 



AM = (M - M C ore), (7) 

where M core is the stellar core mass (j8j* +|8^). If the maximum mass of a star (usually 
it is initial mass of a star, but mass transfer in binary system is able to increase its mass 
over the initial value) M max > 15Af the mass of the core in the main sequence stage 
is determined using l|8j*), and in giant and in supergiant stages using (j8j3). In the Wolf- 
Rayet star stage, if Mwr < 2.5M Q and M max < 2OM it is described using (j8jy), if 
M W r > 2.5M and M max < 2OM as ©), if M max > 2OM using S). 

l-62M°-83, ( a ) 

051+4.21 lg A/ max -0.93(lgJ\/ max ) 2 

M„,„ = { 0.83M& 3 if, (7) (8) 

1.3 + 0.65(Mvra-2.4), (5) 
3.03M^f, (e) 

Scenario C has high mass loss during the WR stage, it may reach 50% of a star 
mass or more here. Mass loss in other stages (MS, giant, supergiant) for stars with 
masses higher than 15M (for less massive stars this scenario is equivalent to a type A 
wind) may reach w 30% of the mass of a star. Total mass loss AM during all stages 
always is larger than in scenario A. 



2.4 Common envelope stage efficiency 

An effective spiral-in of the binary components occurs during the common envelope 
(CE) stage. The effectiveness of the CE-stage is described by the parameter acE = 
AE / AE oro , where AEb = E grav — E t hermai is the binding energy of the ejected 
envelope matter and AE or b is the drop in the orbital energy of the system during spiral- 
in J98). 

acE ( GM a M c _ GM a M d \ = GM d (M d - M c ) 
\ 2a/ 2a,i J Rd 
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where M c is the mass of the core of the mass-losing star of initial mass Md and radius 
Rd (which is simply a function of the initial separation <Zj and the initial mass ratio 
M a /Md, where M a is the mass of the accreting star). 

2.5 Restrictions on key parameters of binary evolution 

Previous estimations of possible value area of parameters of binary stars evolution were 
made in papers ll54ll55l . But during these years some new results have appeared and 
we made some additional calculations. 

Newest observational estimations of neutron star kick during supernova explosion 
is 1 28 1. The authors concluded that characteristic kick velocity vq — 265 km s _1 . 

II II suggested to make more correct estimations for common envelope stage ef- 
ficiency. They argued to make correction of gravitational energy taking into account 
the fact that matter of the stars concentrates to their center GM d(M d -M c ) ^ g ut j n 

Hd A 

calculations they supposed common envelope stage efficiency /iqe to be equal to 1. In 
general, we don't know this parameter exactly. Our coefficient olce is the production 
of the common envelope efficiency [Ice and the parameter A. So we use estimations 
of acE suggested by l54l . 

We would like to note one important circumstance. Poorly known evolution param- 
eters, such as vq, a q , stellar wind, etc. are intrinsic parameters of population synthesis. 
In future they can be defined more exactly or their physical context can be changed: 
distribution of kick velocity might not be Maxwell, complicated hydrodynamics of 
common envelope probably can not be described in terms of a ce and A, distribution 
f(q) might not be power law. So, we have only one way to move ahead in investiga- 
tions, what is to compare predictions of our models with observational data. 

We suggest here two tests for our model: comparisons between calculated number 
of systems of Cyg X-3 type and ratio NN ^ Psr (number of radio pulsar in binary sys- 
tems with a neutron star divided by total number of radio pulsars, binary and single). 
To avoid differences between young and recycled radio pulsar we take only young pul- 
sars. Note that current observational estimation of value of this ratio is ~ 0.001 (more 
exactly, catalogue duplicity rate, |3 1): there are two known young binary radio pulsars 
with a neutron star companion (J2305+4707 and J0737-3039) and > 1500 known ra- 
dio pulsars. As Cyg X-3 type system we take here a black hole with WR-star, mass of 
WR-star is > 7M and orbital period of the pair is less than 10 hours. 

In the Figure 1 we show how ratio NN ^ Fsr (here Nns+ Psr is the calculated num- 
ber of binary neutron stars with radio pulsars and Np sr is the calculated number of 
all radio pulsars) depends on two parameters: kick velocity vq and common envelope 
stage efficiency acE- "Width" of the filled area depicts models with acE in the range 
between 0.2 and 1.0. We used stellar wind type A, a q = for these calculations. 

In the Figure 2 we show OCCO criterion |53] for ratio NN ^^ sr ■ % + §, where O 
is the observed number of the quantity and C is the calculated number of the quantity. 
If O = C this value is equal to 2. As one can see from the Figure 2 high kick velocity 
vq > 200 km s _1 leads to lack of binary radio pulsars with neutron stars. 

In the Figure 3 we present estimated with our model number of Cyg X-3 systems. 
As we can see from this figure, they are practically exclude values of acE < 0.3. 
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Figure 1: This figure shows how ratio NN ^ Fsr (here N^s+Psr is the calculated 
number of binary neutron stars with radio pulsars and Np sr is the calculated number 
of all radio pulsars) depends on two parameters: kick velocity vq and common envelope 
stage efficiency acE- "Width" of the filled area depicts various values of olce hi the 
range between 0.2 and 1.0. 



O/C+C/O, 0=0.0014 
50 i 




Figure 2: This figure shows OCCO criterion 1531 for ratio NN ^ Fsr ■ vq is NS's char- 
acteristic kick velocity. "Width" of the filled area depicts various values of acE in the 
range between 0.2 and 1.0. Note that observation value of this ratio is ~ 0.001. 
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3 Results and discussions 

Frequencies of events which can produce a GRB calculated using "scenario machine" 
are presented the tables [0 and[5] All shown frequencies are normalized to a galaxy 
with mass and star formation rate equal to the mass and the current star formation rate 
in the Milky Way. This suggestion is valid, because even in case of mergers most of the 
events (merging and collapsing) must happen during first billion years 11921 1511 from 
formation of the appropriate systems. It is important to note that the minimal mass of 
the star which evolution remnant is the black hole assumed in the present work to be 
equal to 25 masses of the Sun 1 89 1. 

In the table ^ we have presented frequencies of events which are able to produce 
a GRB. In the first place we put merging of oxygen-neon white dwarfs with oxygen- 
neon, carbon-oxygen and helium white dwarfs and accretion induced collapses (AIC) 
of oxygen-neon white dwarfs in systems with optical companion. Also we showed in 
the tabled frequencies of merging neutron stars with neutron stars and neutron stars 
with black holes. So, in the table^we have showed frequencies calculated using the 
next scenario parameters: stellar wind A, a ce — 0.5, a q — 0, characteristic kick ve- 
locity vq = for neutron stars and black holes. In the table ^frequencies of rotational 
collapses of Wolf-Rayet stars in close binary systems are presented. Since the critical 
period of rotation of WR-star is not reliably fixed yet we have made our calculations 
using two values of critical orbital period P cr u, where P cr n is the maximum orbital 
period of a close binary system in which GRB is able to be formed. Calculations of 
rotational collapses of WR-stars we take into account three types of close binary sys- 
tems: consisting of a black hole and a Wolf-Rayet star (BH+WR), WR-star and main 
sequence star (WR+MS), WR-star and non-degenerated star in Roche lobe overflow 
stage (WR+Rlo). It is necessary to note that in case of stellar wind C minimal period 



10 



Table 1: Frequencies of events in the Galaxy which are able to produce a GRB. Stellar 
wind A, a ce = 0.5, a q = 0. 

Event Frequency, yr _1 

White dwarfs 



ONe+CO 


1.8 


• 10 


-3 


ONe+He 


1.7 


• 10 


-5 


ONe+ONe 


4.9 


• 10 


-4 


ONe AIC 


8- 


io- 


5 


Hyper-nova model 


Pcrit 


= 1 day 




BH+WR 


3.7 


• 10 


-a 


WR+MS 


3.4 


• 10 


-5 


WR+Rlo 


1.1 


• 10 


-5 


Pcrit 


= 3 days 




BH+WR 


9- 


10" 


6 


WR+MS 


2.6 


• 10 


-4 


WR+Rlo 


2.9 


• 10 


-5 



Neutron stars with neutron stars 
and neutron stars with black holes 

mergers, kick velocity vq — 
NS+NS 2.1 • 10" 4 

NS+BH 7.4- 10~ 5 

Comment: P cr it is the minimal orbital period of close bi- 
nary system in which rotational collapse of Wolf-Rayet star 
is possible to produce GRB. 

of a binary star including WR-star just before collapse becomes to about five days; it 
is much higher than estimation of the maximum period of binary in which GRB is able 
to be formed. 

In the tables |2] and 13 one can find how different scenario parameters influence on 
the calculated frequencies of events which we are studying in this work. 

In the table [2] we have presented frequencies of merging of oxygen-neon white 
dwarf with oxygen-neon, carbon-oxygen and helium white dwarfs and accretion in- 
duced collapses (AIC) of oxygen-neon white dwarfs in systems including optical com- 
panion for acE = 1-0. Note that in all cases we calculate low mass (mi < 10M©) 
systems using stellar wind A and a q = 0. Also in the table[2]we have presented fre- 
quencies of rotational collapses of Wolf-Rayet stars in close binary systems which were 
calculated using different scenario parameters. 

In the table[5]one can find frequencies of merging of neutron stars with neutron stars 
and neutron stars with black holes if we used various scenario parameters described 
above. 

During analysis of the frequencies showed in the tables ^ 12 and |3] it is necessary 
to take into account collimation of gamma ray emission in the small space angle |72|. 
It means that only a small part of the GRBs might be observed on the Earth. This fact 
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Table 2: Frequencies of collapses of Wolf-Rayet stars and ONe white dwarfs in the 
Galaxy calculated using different scenario parameters. Stellar wind A. 



Event Frequency, yr 1 


White dwarfs, acE = 


: 1 


ONe+CO 1.1-10 


-3 


ONe+He 7.2 ■ 10 


-5 


ONe+ONe 4.5 ■ 10 


-4 


ONe AIC 2.7 ■ 10 


-5 


Hyper-nova model 


a q = 2, a ce = 0.5 


P crlt = 1 day 


BH+WR 3.2 • 10 


-6 


WR+MS 4.6 ■ 10 


-6 


WR+Rlo 2.2 ■ 10 


-6 


Pcrit = 3 days 


BH+WR 1.1 • 10 


-5 


WR+MS 8.2 ■ 10 


-5 


WR+Rlo 1.1-10 


-5 


a q = 0, a ce = 1 


P crit = 1 day 


BH+WR 2 • 10- 


5 


WR+MS 4.1 • 10 


-5 


WR+Rlo 8.3 ■ 10 


-6 


Peru = 3 days 


BH+WR 3.6 • 10 


-5 


WR+MS 2.9 • 10 


-4 


WR+Rlo 2 • 10- 


5 



Comment: P cr it - the same as in TableQ 
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Table 3: Frequencies of merging of neutron stars with neutron stars and neutron stars 
with black holes under different evolutionary scenario parameters. 



NS+NS, yr- 1 


NS+BH, yr" 1 


Vq, km s 1 


Stellar wind A, a q = 





7.3 ■ 1CT 5 


6.5 • 1CT 5 


50 


2.3 • 1CT 5 


3.7- 10- 5 


100 


7.3 • 1CT 6 


1.7- icr 5 


200 


2.6 • 10" 6 


9.8- 1CT 6 


300 


Stellar wind A, a q = 


2 


4- 10~ 4 


6.1 • 10- b 





1.4- 1CT 4 


5.2 ■ icr 5 


50 


4.7- 1CT 5 


2.8- lO" 5 


100 


1.7 • 1CT 5 


9.4- 10- 6 


200 


6- 1CT 6 


5 • icr 6 


300 


v o = 0, a q = 


Stellar wind C 


NS+NS, yr" 1 


NS+BH, yr" 1 


OtCE 


2 • 1CT 4 


2.3- 10- b 


0.5 


Stellar wind A 


2.4- 1CT 4 


3.9- lQ- b 


1.0 



allows to estimate GRB frequency in the Galaxy: ~ 10 -4 — lO -5 per year 1581 1721 . 
And it is quite possible that collimation depends on the scenario, so table quantities 
can not be directly related to observed frequencies. Comparison of this quantity with 
theoretical estimations of events probably forming GRB flashes (see tabled} leads us 
to a conclusion about potential adequacy of these frequencies for the explanation of 
observable frequency of GRBs. However note that all these frequencies remain rather 
uncertain at present time. Analyzing table |2] it is important to notice considerable 
frequency of ONe degenerated dwarfs merging. For hyper-nova (WR-star collapse) it 
is necessary fiat to use (a q = 0) initial mass ratio distribution and possible maximum 
orbital period limit ~ 1 day. Influence of a model of stellar wind on the frequency 
of NS+NS or NS+BH merging is inessential, but neutron star's or black hole's "kick" 
with magnitude higher than ~ 100 km sec -1 has to be declined if observed GRBs are 
the products of the phenomena studied in this article |48 1. 

We would like also to note that frequencies of WR-stars collapses have no signifi- 
cant differences in cases of critical orbital periods P cr it = 1 day and P cr u — 3 days, 
but if we assume P cr it < 0.5 day all scenarios of GRBs with Wolf-Rayet stars must be 
declined because of zero frequency of such events in frames of our models. 
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